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ABSTRACT 

<\ 

' We investigate analytically a large-scale coherence in the orientation of galax- 

^ I ies embedded in two-dimensional sheet-like structures in the frame of the tidal 

^ ' torque theory. Assuming that the galaxy spin and the surrounding matter fields 

^ ' are intrinsically aligned in accordance with the tidal torque model, we first de- 

^ I rive analytically the probabihty distribution of the galaxy position angles, and 

00 ' evaluate the degree of their inclinations relative to the plane of the sheet. Then, 

^ ' we apply our analytic approach to the nearby spirals in the Local Super Cluster, 

O ■ and provide theoretical explanations about why and to what degree the nearby 

■ spirals are inclined relative to the supergalatic plane. Finally, we conclude that 

I ' the observed large-scale coherence in the orientation of nearby spirals relative 

^ ■ to the supergalactic plane can be quantitatively understood in terms of galaxy 

^ ■ intrinsic alignment predicted by the tidal torque theory, and that the spins of 

• luminous galaxies might be more strongly aligned with the surrounding matter 

r~j ■ than the underlying dark halos. If applied to large scale surveys like Sloan Digital 

j_j ■ Sky Survey (SDSS), our analytic approach will allow us to measure accurately 

■ — ' the strength of galaxy intrinsic alignment which plays a role of statistical error in 
weak lensing searches and can be used as a fossil record to reconstruct cosmology. 

Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

A flurry of recent cosmological research has been focused on studying the galaxy intrin- 
sic alignment. A question of galaxy intrinsic alignment is closely related to the origin of the 
galaxy angular momentum. The standard model for the origin of the galaxy angular momen- 
tum is the tidal torque theory (Peebles 1969; Doroshkevich 1970; White 1984). According to 
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this theory, the galaxy angular momentum is generated by the early tidal interaction with 
the surrounding matter that continued until the moment of the turn-around. 

A generic prediction of the tidal torque theory is the existence of a local correlation 
between the spin and the matter fields: The spin axes of proto-galaxies are preferentially 
aligned with the intermediate principal axes of the tidal shears from the surrounding matter. 

The degree of this galaxy intrinsic alignment must be highest at the moment of the turn- 
around. In the subsequent evolution, the galaxies are likely to lose gradually the initial 
memory about the surrounding matter due to the complicated nonlinear effects, with their 
intrinsic alignments diminishing in strength as a consequence. 

An important question in cosmology is whether the intrinsic alignment remains at 
present epoch to a statistically significant degree or not. It is important in cosmology since 
the galaxy intrinsic alignment plays a role of systematic error in weak gravitational lensing 
searches (Heavens et al. 2000; Croft & Metzler 2000; Catelan et al 2000; Crittenden et al. 
2002; Hirata et al. 2004). Besides, it can provide a crucial clue to the unsolved problem of 
galaxy formation (Dubinski 1992; Catelan & Theuns 1996; Porciani et al. 2002a,b; Navarro 
et al. 2004), and be also used as a fossil record to reconstruct cosmology (Lee & Pen 2000, 
2001, 2002). 

Without knowing the tidal shear field accurately, it is difficult to measure the strength 
of the intrinsic alignment in practice directly from observational data. One possible way 
is to measure it indirectly from N-body simulation data (Lee & Pen 2000; Heavens et al. 
2000; Porciani et al. 2002a,b). However, using N-body data, what one can measure is only 
the intrinsic alignment of underlying dark halos rather than that of the observable luminous 
parts of galaxies. Although the standard theory of galaxy formation assumes that the spin 
axes of luminous parts align with those of dark halos (Mo et al. 1998), it was hinted by 
recent gasdynamical simulations (Chen et al. 2003; Navarro et al. 2004) that the luminous 
parts could conserve the initial memory of the intrinsic shear-spin correlation better than 
the dark matter parts, which prevents us from taking the simulation results as real ones. 
Lee & Pen (2002) attempted to measure the strength of intrinsic alignment directly from 
observational data by reconstructing the tidal shear field with the data from the Point Source 
Catalog of Redshift Survey (PSCz). Although they claimed a detection of intrinsic shear- 
spin correlation, their measurement was far from being accurate due to the high statistical 
noise. 

The difficulty in measuring the galaxy intrinsic alignment lies in the fact that it is a 
local correlation between the galaxy spin and the surrounding matter fields effective only 
over a few Mpc in a three dimensional space. However, it is often observed in the universe 
that the galaxies are surrounded by a coherent two-dimensional sheet-like structure which 
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usually extends over tens of Mpc. For such galaxies that are embedded in a sheet, one may 
expect a large-scale coherent orientation (Navarro et al. 2004). It does not mean that the 
existence of cosmic sheets leads logically to the expectation of a large-scale coherence in the 
orientation of the galaxies. The two concepts are separate: the galaxy intrinsic alignment is 
due to the tidal torques from the surrounding matter, while sheets form through gravitational 
clustering as first collapsed objects (Shandarin et al. 1995). Sheets are just one example of 
the large-scale structures that surround and tidally interact with galaxies. The expectation 
of a coherence in the orientation of the galaxies embedded in a sheet comes not from the 
existence of the sheet itself but from the large-scale coherent pattern in the two-dimensional 
matter distribution of the surrounding sheet. Therefore, it may be possible to predict the 
coherent orientation of the galaxy spin axes relative the plane of the sheet in the frame of 
the tidal torque theory. Furthermore, it may be also possible to evaluate quantitatively the 
strength of this intrinsic inclination of galaxies by an analytic method. The goal of this 
Letter is to explore these possibilities. 



2. ANALYTIC EVALUATIONS 

2.1. Key Assumptions 

To evaluate the intrinsic alignment of galaxies embedded in sheets, we assume the 
following: 

1. The galaxy angular momentum, or the galaxy spin vector (L) is originated by the tidal 
shear effect from the surrounding matter at the proto-galactic stage, which continued 
till the moment of the turn-around. The growth of the proto-galaxy angular momentum 
is determined by its inertia tensor (I) and the local tidal shear tensor (T). The two 
tensors are strongly but not perfectly correlated with each other (Lee & Pen 2000; 
Porciani et al. 2002a,b). The shght misahgnment between the principal axes of the 
two tensors results in a preferential alignment of the proto-galaxy spin axis with the 
intermediate principal axis of the tidal shear tensor (or similarly, the intermediate 
principal axis of the inertia tensor) . 

2. The strength of the galaxy intrinsic alignment with local shears at present epoch can 
be expressed as the following simple quadratic relation (Lee & Pen 2002): 

1 + c 

{LiLj) = — - — Sij — cTikTkj, (1) 

where T is the rescaled traceless shear tensor defined as = Tij/\T\ with = 
Tij — Tr(T)5ij/3, and c is a correlation parameter (see Appendix A in Lee & Pen 2002) 
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introduced to quantify the strength of the intrinsic shear-spin ahgnment in the range of 
[0, 1] . The value of c cannot be determined from 1st physical principles since it includes 
all the nonlinear effects after the moment of the turn-around. One should determine 
its value empirically. 

3. The probability distribution of galaxy spin vectors, -P(L), is well approximated as a 
Gaussian (Catelan & Theuns 1996). 

4. A cosmic sheet is the first collapsed object (Shandarin et al. 1995), forming as early as 
the turn-around moment (Pauls & Melott 1995). The minor axis of T at each point 
on a sheet (i.e., the direction of its maximum compression) is almost perpendicular to 
the plane of the sheet. We approximate the lowest eigenvalue of each T on a sheet as 
zero coherently. 



2.2. Distribution of Position Angles 

According to the third hypothesis in §2.1, the distribution of a galaxy spin vector can 
be written as 



^^^^ " [(27r)3det(M)]V2 



L,(M-i),,L, 



(2) 



with the covariance matrix Mjj = (LiLj) is related to T by equation (1). 



Let Al, A2, A3 represent the three eigenvalues (in a decreasing order) of the tidal shear 
tensor T at some point on the sheet. Following the fourth hypothesis, we approximate 
As's identically as zero at all points on the sheet. Therefore, in the frame of the principal 
axis of the tidal shear tensor at each point, the polar angle ^ of a galaxy spin vector L = 
[L sin 6 cos 0, L sin 6 sin 0, L cos 6) represents the angle of the spin axis relative the plane of 
the sheet (i.e., the galaxy position angle). 

To derive the distribution of the galaxy position angles relative to the plane of the 
sheet, we integrate equation (2) in the principal axis frame of the tidal shear as P{0) — 
Jo Io° -^(-^' ^' 4>)L^dLd4>, and derive the following analytic expression: 



1 A/-, , C2\-^ /^sin2^cos20 sin^^sin^^ cos^ ^ 

P(9) = — (1 + c- 3cA- / ^ H H — d(b 

27rf-LV V Vl + c-3cA? 1 + c - 3cA2 1 + c - 3cA2J 



i=l ^ u \ J- -r ^ - oo/\i J- -r - oL,/\2 -L -r - uo/vg 



1 A/. . „ f9\-5 /"^^ / sin^ ^ cos2 sin^^sin^^ cos^ ^ 

(3) 

where 9 is assumed to be in the range of [0, 7r/2]. Here the three Aj's {i — 1,2,3) are the 
rescaled eigenvalues of T constrained by X^^A? = 1. For the case of A3 = 0, they can be 
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expressed in terms of the two eigenvalues Ai and A2 as 

_ 2Ai — A2 ^ _ 2A2 — Ai _ — (Ai + A2) 

V6(A?-AiA^TA|)' V6(A?-AiA^TAi)' ' " ^6{Xl - AiA^+Af) 

Equations (3)- (4) show how the probabihty distribution of the galaxy position angles, 
P{9), depend on the eigenvalues of Ai and A2 as well as the correlation parameter c. The 
conditional distributions of the two eigenvalues provided that A3 = can be derived by using 

the Baycs Theorem: p(Aj|A3 = 0) = [p(Aj, A3 = 0)/p{\s = 0)] (i = 1,2). Here the two-point 
joint distribution p{\i, A3) and the one-point distribution p{X3) are all obtained analytically 
by Lee & Shandarin (1998). Using these results, we derive 



p(Ai|A3 = 0) = 



p(A2|A3 = 0) = 




where a is the rms fluctuation of the linear density field. 

Figure 1 plots the above conditional distributions (eqs.[5]-[6]). The dotted lines locate 
the most probable values of the two eigenvalues provided that A3 = 0: A™*^^ ~ 0.8a; A™^'^ ~ 
0.3(7. We also determined the conditional mean and the standard deviation of each eigenvalue 
from equations (5)-(6): Ai — 0.86cr, a^^ — 0.28a; A2 = 0.38(7, cr^j = 0.21a. 

Putting A™^ and A™^ into equation (3), we can investigate how the distribution of 
the galaxy position angles relative to the plane of the sheets varies with the value of the 
correlation parameter c. Figure 2 plots P{9) for five different cases of c = 1, 0.7, 0.5, 0.3 and 
(solid, dashed, long-dashed, dot-dashed and dotted lines, respectively). As expected, for 
the case of c = 0, P{9) is uniform, and the larger the value of c is, the more sharply does 
P{9) increase in the small-angle section. In other words, the stronger the intrinsic shear-spin 
correlation is, the more inclined are the galaxy spins relative to the sheet. Equation (3) also 
allows us to evaluate the average position angles: 6 — J^^'^ 6P{9)d6. We find ^ = 29 in unit 
of degree, if c = 0.9 (see in §3). 
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3. THEORY VS. OBSERVATION 

A vivid example of such galaxies that are embedded in a two-dimensional sheet comes 
from the neighborhood of our Milky Way where nearby spirals are assembled in the Local 
Super Cluster (de Vaucouleurs et al. 1991). Although there were some reports of detecting 
a large-scale coherent orientation of the observed nearby spirals relative to the Local Super 
Cluster (e.g., Gregory et al. 1981; Helou & Salpeter 1982; Fhn & Godlowski 1986; Garrido et 
al. 1993), these previous reports often suffer from small sample sizes, hardly being established 
as compelling evidences. 

However, very recently, Navarro et al. (2004) used a relatively large sample of nearby 
galaxies from the Principal Galaxy Catalog (PGC, Paturel et al. 1997), and estimated 
the number distribution of nearby disk galaxies as a function of the supergalatic position 
angles. They found that the spin axes of nearby edge-on spirals {cz < 1200 km/s) are 
strongly inclined relative to the supergalactic plane, and concluded that the observed result is 
consistent qualitatively with the tidal torque theory. Yet, they did not make any quantitative 
comparison of the observed result with the predictions of the tidal torque theory. 

To make a quantitative comparison, we evaluate the number distribution of galaxies 
embedded in sheets with equation (3). Figure 3 plots the results. The (red) dot-dashed line 
represents our theoretical prediction with the choice of c = 0.9 (the best-fit value), while the 
histogram corresponds to the observational result from PGC (see Fig. 2 in Navarro et al. 
2004). Obviously, the theoretical curve agrees with the histogram very well. Navarro et al. 
(2004) determined the average value of the position angles from the observed edge-on disk 
galaxies to be Bobs ~ 25 in unit of degree, which is in agreement of the theoretical value that 
we obtained assuming c = 0.9 in §2. 

It is worth noting that the best-fit value of c = 0.9 is three times higher than that found 
in N-body simulation (Lee & Pen 2000, 2002). It may imply that the luminous parts of 
galaxies tend to keep the initial memory of the surrounding matter field better than their 
dark matter counter parts, consistent with the results from recent gasdynamical simulations 
(Chen et al. 2003; Navarro et al. 2004). 

4. CONCLUSION 

We have predicted analytically a large-scale coherence in the orientation of galaxies em- 
bedded in sheets using the tidal torque theory. Our analytic model reproduces the observed 
inclinations of nearby spiral galaxies in the Local Super Cluster remarkably well, providing 
a physical quantitative understanding of the observables. It should be, however, emphasized 
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here that the prediction of the coherent orientation of galaxies embedded in a sheet should 
not be taken as identical to the existence of the sheet itself. The coherent orientation of the 
galaxies is caused by the tidal torquing from the surrounding matter, which is different from 
the formation of cosmic sheets in the universe. 

By comparing our analytic prediction with the observational result, we also measured 
quantitatively the strength of the galaxy intrinsic alignment, and found that the luminous 
galaxies seem to conserve their angular momentum better than the underlying dark halos 
after the moment of the turn-around. Our final conclusion is this: our analytic calculation, if 
applied to observational data from large scale surveys like Sloan Digital Sky Survey (SDSS), 
will allow us to measure accurately the strength of galaxy intrinsic alignment that will have 
an impact in weak gravitational searches, and can also be used as a fossil record to reconstruct 
cosmology. 

We thank C. Park and U. Pen for useful discussions. We also thank the referee for 
helpful suggestions. This work is supported by the research grant of the Korea Institute for 
Advanced Study (KIAS). 
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Fig. 1. — Conditional probability density distributions of the largest and the second largest 
eigenvalues of an intrinsic shear tensor provided that its smallest eigenvalue has the value of 
zero. The dotted line locates the most probable eigenvalue for each case. 
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Fig. 2. — Probability density distributions of the angles of the minor axes of galaxies em- 
bedded in a two-dimensional sheet relative the plane of the sheet. The solid, dahsed, long- 
dashed, the dot-dashed, and dotted lines represent the cases of the correlation parameter 
c — 1, 0.7, 0.5, 0.3 and 0, respectively. 
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Fig. 3. — Number distributions of supergalactic position angles of spirals embedded in a 
sheet. The dot-dahsed curve represents the theoretical prediction assuming c = 0.9, while 
the histogram is the observational result from PGC catalog (see Fig. 2 in Navarro et al. 
(2004) 



